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CHAPTER  I 


INTRODUCTION 

Radar  Fundamental a 

In  it*  siaplaat  fora,  a  radar  ayataa  transmit*  a  signal 
and  wait*  for  a  raflactad  aignal  (acho)  froa  a  poariblt 
targat.  If  an  acho  doa*  exist,  it  will  ba  aiailar  to  tha 
transmitted  aignal,  though  exhibiting  aoaa  diffarancaa. 
Inforaation  about  tha  targat* a  location,  aiaa,  velocity,  and 
direction  of  aovaaant  can  ba  extracted  froa  theae  diffarancaa. 

If  on  tha  other  hand,  there  i*  no  acho,  it  ia  aaauaad  that  no 
targat  exiata  (in  that  direction)  and  tha  aignal  i*  tranamittad 
in  another  direction  (Txannea  1985). 

Tha  radar  aquation  ia  uaad  to  determine  tha  aaximua  range 
at  which  a  targat  can  ba  detected.  If  tha  ayataa  eaita  a  aignal 
with  initial  average  power  of  P^,  watta  and  aaata  an  object  at 
a  diatanca  R  froa  tha  radar,  tha  return  power  of  tha  acho,  PR, 
i*  given  by 


P 


R 


(1) 


where  K  depend*  on  various  other  parameters  of  the  system  such 
as  target  cross  section  and  antenna  gain.  A  useful  fora  cf  r.he 


2 


radar  aquation  ia  found  by  denoting  tha  ainiaum  detectable  power 
by  S  in,  and  aolving  for  R, 


max 


(sf-> 

min 


(2) 


where  R  rapraaanta  tha  aaxiaua  rang#  of  tha  ayataa  (Tsannaa 
1983) . 

Rang a  information  ia  found  by  tranaaitting  a  pulaa  and 
aaaauring  tha  tiaa  dalay,  At,  between  tha  tranaaittad  and 
racaivad  pulaa.  Sinca  alaetroaagnatic  anargy  travala  at  tha 
apaad  of  light,  c,  tha  diatanca  to  tha  targat  will  ba 


R  - 


cAt 

2 


(3) 


where  tha  factor  of  2  accounta  for  tha  round  trip  (Cook  and 
Barnfald  1967). 

To  measure  thia  dalay,  a  diatinct  point  of  rofaranca  ia 
raquirad.  Though  tha  baginning  or  and  of  tha  pulaa  may  hava 
aharp  adgaa  for  reference,  thaaa  adgaa  bacoaa  roundad  or  obacurad 
whan  band-limited  or  in  tha  praaanca  of  noiaa.  To  minimise  tha 
ranga  error,  tha  tiaia  croaa-corralation  function  of  tha  incoming 
echo  signal  with  tha  original  pulaa  ia  takan  using  a  Batched 
filter,  which  naxinisas  tha  signal  to  noiaa  ratio  (SNR)  in  the 
presence  of  white  Gaussian  noise.  Thia  iapliea  that  a  d-**‘  net 
point  of  reference  oust  ba  present  at  tha  output  of  tha  matched 
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filter.  Coed  rang*  resolution  requires  the  output  of  the  notched 
filter  to  have  a  high  sharp  peak,  and  only  a  single  peek,  so  thet 
a  number  of  targets  can  be  distinguished  fro*  one  another 
(Rihacsek  196$). 

A  very  narrow  pulse  would  possess  the  characteristics  needed 
for  good  range  resolution,  but  due  to  power  limitations  in  the 
transmitter,  the  maximum  range  detection  capability  would  be 
limited.  Bence,  in  order  to  meet  range  requirements,  the  radar 
system  designer  seems  to  be  faced  with  a  trade-off.  This 
trade-off  is  reduced  using  pulse  compression  techniques,  which  is 
another  significant  advantage  of  using  matched  filters.  The 
pulse  can  actually  be  made  as  wide  as  necessary  to  meet  range 
requirements,  then  coded  with  wideband  modulation  information  to 
meet  the  range  resolution  requirement  (Cook  and  Bernfeld  1967). 

If  a  target  is  moving  toward  a  radar  systam,  then 
the  frequency  of  the  return  signal  will  be  higher  than  that  of 
the  transmitted  signal.  Similarly,  if  the  target  is  moviig 
away  from  the  radar,  the  frequency  of  the  return  signal  will  be 
lower  than  that  of  the  transmitted  signal.  This  phenomenon  is 
known  as  the  Doppler  effect  and  the  amount  of  change  in  frequency 
(fg),  is  used  to  compute  the  target's  radial  velocity  (\>r).  If 
is  much  less  than  the  speed  of  light,  the  doppler  frequency  is 
approximated  by 


\fti  if*  Muni)  *  ys  a‘PL»  »mSms  MmxM a  Ml 


*h«M  f t  represents  the  transmission  fraquancy  (Wheeler 
1967). 

Linaar  FW 

The  most  coaaaon  font  of  pulaa  compression  uaad  in  modern 
radar  ayataaia  la  linaar  frequency  modulation  (chirp),  which  la 
raallaad  by  linearly  varying  tha  carrier  fraquancy  of  tha  FN 
waveform  during  tha  Interval  of  tha  pulaa.  Chirp  radar  provldaa 
good  range  and  velocity  measurements,  aa  wall  aa  roaolution 
(Tsannea  1963).  Range  raquiraaanta  (within  limits)  can  be  aat  by 
increaaing  the  width  of  tha  tranaaiittad  pulaa.  Uaing  Batched 
filters,  thia  can  be  accoapliahed  without  sacrificing  bandwidth. 

Though  it  would  seea  that  chirp  offara  all  tha  daaired 
character la tics  of  a  radar  aignal,  thia  la  not  the  caaa.  It 
measures  range  vary  wall  whan  tha  target  la  not  aoving  (no 
Doppler  ahifta),  but  when  tha  target  la  aoving,  range  cannot  bo 
calculated  directly.  Though  thia  range  error  can  be  accounted 
for  when  there  la  only  one  target,  it  becoaea  a  problea  when 
there  are  a  number  of  targeta  within  one  pulaa  period  with 
differing  velocitiea.  Therefore,  chirp  radar  la  aoat  auitable 
for  applicationa  where  the  expected  differences  in  velocity  are 
relatively  saall  during  any  singla  pulaa  period  (Cook  and 
Dernfeld  1967).  Another  disadvantage  of  tha  chirp  radar  is 'that 
the  autocorrelation  tine  aidelobea  ore  relatively  high,  which 
could  cause  range  resolution  problems  (Booher  1985). 
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Non-Linear  FM 

Many  approaches  have  been  investigated  in  an  effort  to 
minimize  the  problems  inherent  to  linear  FM.  One  such  approach 
is  to  vary  the  frequency  of  the  transmitted  pulse  in  a  non-linear 
fashion  (non-lincar  FM) .  Work  done  by  Booher  suggests  there  may 
be  some  merit  in  varying  the  pulse  frequency  using  the  Blackman 
function  as  the  modulating  waveform.  He  showed  that  under 
certain  conditions,  this  new  non-linear  FM  could  perform  better 
than  linear  FM  when  range  reaolution  is  of  primary  concern 
(Booher  1983).  This  pulse  is  further  investigated  to  see  if  a 
wider  range  of  velocities  can  be  detected  without  introducing 
significant  range  error  due  to  Doppler  shifts. 

SAW  Devices 

Surface  acoustic  wave  (SAW)  technology  offers  a  means  of 
processing  complex  waveforms  onto  devices  that  are  much  smaller 
and  more  reliable  than  previous  techniques.  The  planar  nature 
of  SAW  devices  allows  them  to  be  fabricated  using  standard 
lithography  techniques  used  by  the  semiconductor  industry.  This 
process  is  highly  repeatable  and  relatively  inexpensive. 

Impulse  response  model  design  techniques  (Hartmann,  Bell 
and  Rosenfeld  1973)  are  used  in  order  to  provide  a  straight¬ 
forward  appTvach  to  producing  sample  devices  which  implement  buth 
linear  FM  and  the  new  non-linear  FM  pulse  compression  waveforms. 


Design  considerations  and  implementation  procedures  are  presented 
in  detail. 


CHAPTER  II 


OBJECTIVE  OF  PROPOSED  WORK 

The  objective  of  this  thesis  is  to  further  investigate 
the  characteristics  of  the  new  non-linear  FM  waveform  introduced 
by  Booher.  SAW  devices  will  be  designed  and  photomasks  will  be 
generated  using  both  linear  FM  and  the  new  non-linear  FM 
function.  Additionally,  V-chirped  devices  using  both  linear  FM 
and  the  new  non-linear  FM  function  will.be  investigated.  Design 
considerations  will  be  discussed  and  computer  aided  design 
software  will  be  written.  Comparisons  will  be  made  on  the 
impulse  response,  frequency  response,  and  matched  filter 
responses  of  the  filters.  The  Doppler  effect  on  both  linear 
FM  and  the  new  non-linear  FM  waveform  will  be  simulated. 
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CHAPTER  III 


WAVEFORM  DESCRIPTIONS 

Using  impulse  response  design  techniques,  SAW  device  design 
is  straightforward  if  the  time  waveforms  can  be  sampled 
accurately  (Hartmann,  Bell,  and  Rosenfeld  1973).  This  chapter  is 
concerned  with  the  derivation  of  the  general  form  of  the  linear 
and  new  non-linear  FM  waveforms.  Though  up-chirped,  down-chirped 
and  V-chirped  SAW  devices  are  considered,  symmetry  arguments  are 
used  (Chapter  IV)  so  only  the  up-chirped  forms  of  the  equations 
are  necessary.  Additionally,  plots  of  the  spectra,  matched 
filter  outputs,  and  spectra  of  the  matched  filter  outputs  for 
both  linear  FM  and  the  new  non-linear  FK  are  presented.  With  the 
exception  of  figures  1  and  2,  the  plots  provided  in  this  chapter 
were  accomplished  using  the  FFT  and  graphic,  capabilities  of  the 
surface  acoustic  wave  computer  aided  design  program  (SAWCAD) 
developed  at  the  University  of  Central  Florida  (Malocha  and 
Richie  1984). 


FM  Waveforms 

In  general,  an  FM  waveform,  Xp^ft),  can  be  represented  by 
the  equation  (Ziemer  and  Tranter  1985) 


t 

Xp^Ct)  -  A  cos  [wct  +  <Kf  /  mOOdT  +  <j>o)]  (5) 

o 
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where : 

i#(T)  •  the  modulating  signal 
A  ■  the  peak  amplitude  of  carrier 

wc  ■  the  angular  frequency  when  m(t)  is  zero 

m  the  frequency  deviation  constant 
4>o  «  the  initial  phase 

The  phase  deviation,  $(t),  of  XpM(t)  is 

t 

<Kt)  -  K,  /  m(T)dT  +  4> 


(6) 


The  instantaneous  frequency,  o^Ct),  of  XpH(t)  is  found  by  taking 
the  derivative  of  the  bracketed  term  in  equation  (5)  with  respect 
to  time  which  yields 


(<^(0  -  u>c  +  Kf  m(t)  (7) 

Linear  FM 

Linear  FM  is  an  FM  waveform  whose  instantaneous  frequency 
changes  in  a  linear  fashion  with  respect  to  time.  Figure  1  shows 
the  modulating  signal  used  to  accomplish  this.  Note,  the 
remaining  equations  (6-16)  in  this  chapter  are  only  valid  from  0 
to  T  since  they  represent  pulses,  not  continuous  waveforms.  The 
linear  FM  modulating  signal,  m^p^t),  is  given  by  the  equation 

“LFM^  "  T~ 


(8) 
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The  instantaneous  frequency  of  the  linear  FM  waveform,  wlFM< 
is  found  by  substituting  equation  (8)  into  equation  (7)  to  yield 

aK  t 

'  "c  *  —  <» 

The  phase  deviation  of  the  linear  FM  waveform,  pn < t ) ,  is  found 

by  substituting  equation  (8)  into  equation  (6)  and  performing  the 
integration  such  that 

2 

*LFM(t^  "  2T~  t  +  *0  (10) 


Now,  the  linear  FM  waveform,  K^FM'^*  can  be  written  by 
substituting  equation  (10)  into  equation  (5),  which  yields 


W*> 


aKf  2 

A  cos  (wct  +  t  +  <J>q) 


(11) 


Non-Linear  FM 

This  new  non-linear  FM  uses  a  shifted  form  of  the  Blackman 
function  as  the  modulating  signal  of  the  FM  waveform  (Booher 
1985).  Shown  in  Figure  2,  the  modulating  signal  of  the  new  non¬ 
linear  FM  waveform,  mur_u(t),  is  represented  by 

NLr  W 


m^p^t)  -  a{0.43  +  0.5  cos^-'^j  +  0.07  cos[ <12) 
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The  instantaneous  frequency  of  the  new  non-linear  FM  waveform, 

1*  found  by  aubatituting  equation  (12)  into  aquation 
(7)  to  yield 

WNLIM<t)  "  Wc  +  ■V0,43  +  0,5  coal-^^-]  +  0.07  cos [&&=&] )  (13) 

The  phaae  deviation  of  the  new  non-linear  FH  waveform,  FM(t) , 
is  found  by  aubatituting  equation  (12)  into  equation  (6)  such 
that 


<*>NLFM(t) 


(14) 

aKf  /  (0.43  +  0.5  costal]  +  0.07  cos  [241=31,  )(lt  +  «, 

t  Q  1  l  -O 


Performing  the  required  integration  yields 


(15) 


ViM(t)  "  «M°*43t  +  ^r1  8in[27T^-~-^-]>  +  4> 


2tt 


Finally,  the  new  non-linear  FM  waveform,  X^^p^(t) ,  can  be  written 
by  substituting  equation  (15)  into  equation  (5),  which  yields 

XNLFM(t)  “  A  cos|a,ct  +  aKf{°’43t 

+  sin[-^=^-)  +  sin[-^^]}  +  <t>oj 

Design  Characteristic  Waveforms 
In  general,  it  is  very  difficult  (if  not  impossible)  to 
physically  see  the  varying  frequency  of  practical  FM  waveforms. 


(16) 


\»\wuwuw\AruvwiAcuinjw^^nAAJWUwuuuw\  wuAiuukAMAMAKrTHAiiRiLA'iLn’MJiYAMJi-MA  wirtLn  tunutunut  ii*v.TMJ!iijr«jnL<kJtwnu 
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For  this  reason,  computer  plots  of  tho  actual  designed  chirpad 
waveforms  will  not  ba  included,  though  a  hypothetical  exanple  of 
an  up-chirped  linear  FN  wave fonn  (Figure  3)  and  a  V-chirped 
linear  FM  waveform  (Figure  4)  has  been  included  for  pedagogical 
purposes.  The  waveform  of  Figure  3  chirps  from  40  NHs  to  100  MHx 
in  approximately  0.1  uaec,  while  the  waveform  of  Figure  4  chirps 
from  40  NHs  to  100  NHs  then  back  to  40  NHs  in  approximately  0.1 
usee.  The  waveforms  of  the  actual  designs  chirped  from  60  NHs  to 
80  NHs  in  approxinately  3  usee. 

Frequency  Response  of  the  Filters 

The  spectrum  of  the  linear  FN  waveform  can  be  represented  in 
closed  form  (Cook  and  Bernfeld  1967).  Normalised  amplitude  plots 
of  the  linear  and  new  non-linear  FN  spectra  are  shown  in  figures 
3  and  6,  respectively.  These  plots  were  accomplished  by 
performing  en  FFT  on  the  actusl  time  waveforms  used  in  the  SAW 
device  design. 

The  frequency  responses  of  the  receiving  (matched)  filters 
were  found  by  multiplying  the  impulse  response  of  the 
trensmitting  filter  with  that  of  its  corresponding  matched 
filter.  Figures  7  and  8  show  the  frequency  responses  of  the 
matched  filter  outputs.  It  should  be  noted  that  since  the  new 
non-linear  FN  actually  has  two  peaks,  linear  FM  would  perform 
better  in  terms  of  velocity  resolution. 


IMUdUKM(B19UttUtfMIUlAJ0UOU(lti(N(NlfUBUUytAJ(JU(Nl^UUlAAAM  *AA*nii^\K»U<nXT\xx*TVTLj^  IK*  ftkTUflLtfwn) 
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Tiae  Reaponaea  of  tho  Matched  Flltora 
Tho  tlaa  roaponaa  of  tho  matched  filtora  (tho 
autocorrelation  function)  ia  found  by  taking  the  inverae  FFT  of 
the  corraaponding  matched  filter  frequency  roaponaa.  Matched 
filter  outputa  for  linear  and  the  new  non-linear  FM  wavefoma  are 
praaantad  in  figuraa  9  and  10  raapectivaly.  The  aidelobe  level 
of  the  linear  FM  ia  approximately  13,7  dB  down,  while  that  of  the 
non-linear  FM  ia  only  7,6  dB  down.  In  thia  case,  the  linear  FM 
providea  better  range  reaolution. 


CHAPTER  IV 


DOPPLER  SIMULATION 

The  Doppler  effect  cen  be  simulated  on  e  computer  by 
shifting  the  center  frequency  of  the  transmitted  pulse  and 
processing  the  resulting  signal  through  the  original  matched 
filter.  This  was  accomplished  by  taking  the  FFT  of  the  frequency 
shifted  up-chirped  waveform,  and  multiplying  it  with  the  FFT  of 
the  original  down-chirped  waveform.  An  inverse  FFT  war  ".hen 
performed  on  the  product  to  yield  the  desired  Doppler  shifted 
matched  filter  response.  The  simulation  was  performed  on  both 
linear  FM  and  the  new  non-lineaz  FM  waveforms  for  several  Doppler 
frequency  shifts  ( f  g) ,  and  on  the  V-chirped  versions  of  the 
waveforms . 

Distortion  effects  on  the  output  of  the  original  matched 
filter  caused  by  Doppler  shift  conditions  are  loss  in  peak 
amplitude  and  a  time  shift  of  the  waveform.  For  linear  FM,  the 
amplitude  degradation  is  essentially  bounded  by  a  triangle  with 
its  baseline  extending  from  -T  to  +T.  This  triangle  is  the  auto¬ 
correlation  of  the  rectangular  envelope  used  and  this 
relationship  is  a  property  of  the  linear  FM  waveform.  For  fg 
less  than  half  the  total  chirped  bandwidth,  a  good  approximation 
for  the  time  shift,  t  ,  of  the  linear  FM  waveform  is  (Cook  and 

o 

Bernfeld  1967) 
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Both  of  these  distortion  affects  can  be  easily  verified  from  the 
linear  FM  Doppler  simulation  presented  here. 

All  of  the  following  plots  are  of  the  matched  filter  output 
for  the  Doppler  frequency  indicated.  The  center  frequency  of  the 
receiving  filter  was  set  at  70  MHz,  the  dispersion  time  at  3 

usee,  and  the  bandwidth  at  20  MHz.  The  center  frequency  of  the 
Doppler  shifted  pulse  was  set  at  70  MHz  plus  the  indicated 
Doppler  frequency  shift,  the  dispersion  time  at  3  usee,  and  the 
bandwidth  at  20  MHz. 

Figures  11  and  12  are  for  a  Doppler  frequency  shift  of  0, 
for  linear  FM  and  the  new  non-linear  FM  respectively.  The 
following  plots  in  this  simulation  are  scaled  to  the  peak 
amplitudes  of  these  two  in  order  to  show  the  relative  amplitude 
degradation  dee  to  the  Doppler  shifts. 

Figures  13,  15,  17,  19,  and  21  represent  Doppler  frequency 
shifts  of  +1,  +3,  +5,  +7,  and  -5  MHz,  respectively,  for  the 
linear  FM.  The  resulting  Doppler  distorted  waveforms  are  as 
predicted  by  the  triangular  envelope  for  amplitude  degradation 
and  equation  (17)  for  the  time  shift.  The  results  are  summarized 
in  Table  1. 

The  Doppler  distorted  waveforms  for  the  new  non-linear  FM 
are  presented  in  figures  14,  16,  18,  20,  and  22  for  the  same 
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Figure  16.  New  non-linear  FM  matched  filter  output  with 
Doppler  frequency  shift  of  3  MHz. 
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Figure  17.  Linear  FM  Matched  filter  output  with 
Doppler  frequency  shift  of  5  MHz. 
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Figure  18.  New  non-linear  FM  matched  filter  output  with 
Doppler  frequency  shift  of  5  MHz. 
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Figure  19.  Linear  FM  matched  filter  output  with 
Doppler  frequency  shift  of  7  MHz. 
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Figure  20.  New  non-linear  FM  matched  filtor  output  with 
Doppler  frequency  shift  of  7  MHz. 
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Figure  21.  Linear  FM  matched  filter  output  with 
Doppler  frequency  shift  of  -5  MHz. 
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figure  22.  New  non-linear  FM  matched  filter  output  with 
Doppler  frequency  shift  of  -5  MHz. 


respective  Doppler  frequency  shifts  used  In  the  lineer  FH 
simulation.  These  results  ere  summarised  in  Tsble  2. 
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TABLE  1 

DOPPLER  DISTORTION  EFFECTS  ON  LINEAR  FM 


DOPPLER  FREQUENCY 
(MBs) 

AMPLITUDE  DEGRADATION 
<dB) 

TIME  SHIFT 
(usee) 

1 

0.45 

-0.15 

3 

1.41 

-0.45 

5 

2.30 

-0.75 

7 

3.74 

-1.05 

-5 

2.50 

+0.75 

TABLE  2 

DOPPLER  DI STATION  EFFECTS  ON  THE  NEW  NON-LINEAR  FM 


DOPPLER  FREQUENCY 
(MHs) 

AMPLITUDE  DEGRADATION 
(dB> 

TIME  SHIFT 
( usee ) 

1 

3.84 

-0.12 

3 

8.12 

-0.30 

5 

9.39 

-0.48 

7 

9.86 

-0.68 

-5 

9.39 

+0.48 

Though  the  amplitude  degradation  due  to  Doppler  shifts  is 
much  more  severe  for  the  new  non-linear  FM,  a  significant 
improvement  is  seen  in  terms  of  time  shift  distortion.  This 
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characteristic  could  be  exploited  to  yield  e  systea  which  ia 
capable  of  detecting  e  larger  range  of  velocitiea  at  the  expenae 
of  more  power  in  the  tranaaitted  pulae. 

Another  aignificant  characteristic  ia  that  the  aidelobea  on 
one  aide  (the  aide  toward  the  aero  velocity  point)  for  the  new 
non-linear  FM  wavefora  have  virtually  diaappeared  when  aubjected 
to  Doppler  ahifta.  This  provider  information  not  iaaediately 
apparent  in  the  linear  FM  wavefora.  Firat  of  all,  it  aervea  ea  a 
flag  to  indicate  whether  a  target  ia  aoving  or  not.  Secondly,  if 
the  target  ia  aoving,  it  indicatea  whether  it  ia  aoving  toward  or 
away  froa  the  radar. 

It  ahould  be  noted  that  equation  (17)  appliea  to  a  ayataa 
in  which  the  tranaaitted  pulae  waa  an  up-chirp  and  the  aatched 
filter  iapulae  reaponae  waa  a  down-chirp.  If  the  ayatea  were 
reveraed  such  that  the  tranaaitted  pulae  waa  down-chirped,  the 
tine  shift  would  be  in  the  opposite  direction.  If  the 
tranaaitted  pulse  chirped  up  and  down  (V-chirp),  then  there  would 
be  a  shift  in  both  directions. 

The  V-chirp  offers  a  Beans  of  aeasuring  target  velocity 
without  the  need  for  a  bank  of  narrow  band  filters.  When  a 
V-chirped  pulae  aeets  a  aoving  target,  the  resulting  matched 
filter  output  has  two  naxiaa  due  to  the  Doppler  effect  which  are 
aymaetric  about  the  aero  velocity  point.  This  essentially 
creates  its  own  reference  point  from  which  the  time  shift  can  be 


29 


aaaaurad  and  tha  velocity  coaputad.  Actually,  tha  V-chirp 
aaaauraa  spaad,  not  valocity,  ainca  thara  la  no  way  of  knowing 
for  a  givan  aaxiaua,  whathar  it  waa  eauaad  by  tha  up-  or  down- 
chi  r  pad  portion  of  tha  pulaa,  thua  loaing  tha  aanaa  of  targat 
diraction.  Alao,  dua  to  tha  praaanea  of  two  aaxiaa,  targat 
rasolving  capability  la  liaitad. 

Figuraa  23  and  24  ara  for  a  Dopplar  fraquancy  ahift  of  0, 
for  linaar  FN  and  tha  naw  non-linaar  FN  raapactivaly,  Figuraa 
25  and  26  rapraaant  a  Dopplar  fraquancy  ahift  of  -5  NHs  and  ara 
acalad  to  tha  paak  aaiplituda  of  thair  raapactiva  0  Dopplar 
fraquancy  ahift  wavaforaa  in  ordar  to  ahow  tha  ralativa  aaiplituda 
degradation.  Table a  3  and  4  auaaarisa  tha  raaulta  of  tha  Dopplar 
aiaiulation  for  tha  linaar  V-chirp  and  tha  naw  non-linaar  V-chirp, 
raapactivaly . 


TABLE  3 

DOPPLER  DISTORTION  EFFECTS  ON  V-CHIRPED  LINEAR  FN 


DOPPLER  FREQUENCY 

AMPLITUDE  DEGRADATION 

TINE  SHIFT 

(NHs) 

(dB) 

(uaac) 

0 

0.00 

0.00 

-5 

8.92 

0.40 
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Figure  23.  V-chirped  linear  FN  Matched  filter  output 
with  no  Doppler  frequency  ahift. 


Figure  24.  V-chirped  non-linear  FN  matched  filter  output 
with  no  Doppler  frequency  ahift. 
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Figure  25.  V-chirped  linear  FN  notched  filter  output 
with  Doppler  frequency  shift  of  -5  MHt. 
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Figure  26.  V-chirped  non-lineer  FM  notched  filter  output 
with  Doppler  frequency  shift  of  -5  MBs. 
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TABLE  A 

DOPPLER  DISTORTION  EFFECTS  ON  V-CHIRPED  NON-LINEAR  FN 


DOPPLER  FREQUENCY 

AMPLITUDE  DEGRADATION 

TIME  SHIFT 

(KBs) 

(dl) 

(usee) 

0 

0.00 

0.00 

-3 

11.66 

0.26 

As  before,  linear  FN  off art  laaa  amplitude  degradation, 
while  the  new  non-linear  FN  provide*  laaa  tine  shift  distortion. 
For  V-chirp,  this  tine  shift  distortion  is  actually  desired, 
since  it  is  used  to  conpute  velocity.  Also,  no  range  error  is 
introduced  when  the  ■ id point  between  the  two  naxiaa  is  taken  aa 
the  reference.  Again,  if  the  need  for  a  nora  dynanic  range  of 
velocities  exists,  perhaps  the  new  non-linear  FN  should  be 


considered 


CHAPTER  V 


SAW  DEVICE  DESIGN  CONSIDERATIONS 

The  linear  FN  and  the  new  non-linear  FM  tine  waveforms  can 
be  completely  deacribed  mathematically  in  a  relatively  simple 
form.  In  contrast,  the  mathematical  representation  of  the 
frequency  spectrum  of  these  signals  is  very  complex  and 
approximations  must  be  used  in  the  derivation.  Impulse  response 
model  design  of  SAW  devices  lends  itself  very  well  to  this 
situation  because  of  the  correspondence  between  the  location  of 
the  electrodes  on  the  transducer  and  the  signal  generated  by  an 
impulse  of  acoustic  energy  traveling  under  the  electrodes 
(Hartmann,  Bell,  and  Rosenfeld  1973). 

Critical  Points  of  the  Chirped  Waveform 

The  proper  positioning  of  the  transducer  electrodes  is 
dependent  on  the  location  of  the  positive  peaks  (peaks), 
negative  peaks  (valleys),  and  nulls  of  the  waveform.  The 
distance  between  electrodes  will  vary  proportional  to  the 
changing  frequency. 

Closed  Form  Solution 

For  the  linear  FN  waveform,  the  critical  points  can  be  found 
in  closed  form.  By  normalizing  equation  (11)  and  allowing  the 


initial  phase,  <Pq ,  to  ba  -tt/2,  tha  following  form  of  the 
aquation  is  realised 


yi  -  •*>=  <v +  f t2) 


(18) 


where  V  -  aKf/T  is  the  chirp  slope.  The  peaks  of  the  waveform 
are  found  when  X^p^(t)  i*  equal  to  1.  Allowing  t^  to  represent 
the  time  at  which  tha  peaks  occur  the  equation  becomes 


sin  (u>  t  +  y  t  2)  -  1 

C  p  l  D 


(19) 


which  can  be  written  as 

2  tp2  +  Vp  "  8in”1^1>  "  0  (20) 

or  simply 

2  fcp2  +  Vp  ”  +  2n7T)  "  0  (21> 

where  n  is  an  integer  and  accounts  for  multiple  peaks.  Using  the 
quadratic  formula  to  solve  for  t^  and  ignoring  the  negative  root 
yields 


J  2 

-  0J_  +  o)_  +  uir(4n  +  1) 


35 


The  valleys  and  nulls  of  the  waveform  can  be  found  in  a  aim  Liar 
matter  by  setting  X^pg(t)  in  equation  (18)  equal  to  -1  and  0, 
respectively.  The  results  of  doing  so  gives 


t 

v 


-  a)  +  4  *  +  vm(4n  “  1) 
c  c 

y 


(23) 


and 

/  2 

-  ui  +  u>„  +  2ymr 

t  _  - £_ - £ -  (24) 

n  y 

where  ty  and  tR  represent  when  the  valleys  and  nulls  of  the 
waveform  occur. 

A  closed  form  solution  for  the  critical  points  of  the  new 
non-linear  FM  waveform  is  not  as  easy.  From  equation  (16),  it 
can  be  seen  that  finding  the  critical  points  is  no  longer  a 
matter  of  solving  a  simple  quadratic.  It  was  suggested  by  Booher 
(1985)  that  this  task  could  be  accomplished  using  an  iterative 
approach  of  oversampling.  Although  this  approach  is  valid,  a 
more  efficient  way  of  doing  this  is  possible. 

The  Walkerdid  Algorithm 

Although  this  algorithm  was  developed  for  finding  the 
critical  points  of  the  new  non-linear  waveform  discussed  in  this 
thesis,  it  can  also  be  applied  to  linaar  FM  or  other  non-linear 
FM  waveforms.  The  algorithm  is  applied  to  the  up-chirped 
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(increasing  frequency)  fora  of  the  FM  equation.  Symmetry  is  used 
if  the  desired  critical  points  were  for  tho  down-chirped  or 
V-chirped  fora  of  the  equation. 

First,  define  Td  as  the  desired  amount  of  dispersion  time  in 
the  waveform  and  Ta  as  the  actual  amount  of  dispersion  time  it 
takes  t  >  go  from  the  minimum  to  maximum  frequency  and  such  that 
the  final  waveform  starts  on  a  null  and  has  positive  slope  and 
ends  on  a  null  and  has  negative  slope  as  shown  in  Figure  27.  For 
up-  or  down-chirped  waveforms,  Ta  is  approximately  equal  to  Td 
and  for  V-chirped  waveforms,  Ta  is  approximately  equal  to  Td/2. 
Letting  Xu<t)  represent  the  up-chirped  waveform,  X^(t) 
represent  the  down-chirped  waveform,  and  Xy(t)  represent  the 
V-chirped  waveform,  symmetry  suggests 

X.(t)  =  X  (Ta  -•  t)  (25) 

a  u 

end 

!X(t)  0  <  t  <  Ta 

(26) 

Xu(2Ta-t)  Ta  <.  t  <.  2Ta 

For  up  or  down-chirped  linear  FM,  the  value  of  Ta,  (Ta  ), 

LFM 

is  found  by  evaluating  equation  (16)  for  t  ■  Ta  and  setting 

LFM 

it  equal  to  sere  (want  to  end  on  a  null)  such  that 


.in  U«c  +  f  *  0 


(27) 
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Now,  noting  that  V  *  aK^/T  ■  aK^/Ta^.^,  •^uation  (27)  is  written 

aa 

•Kf 

sin  [<wc  +  ~ r)  T^m]  -  0  (28) 


r 

or 


•Kf 

TaLEW(wc  +  ~ )  “  nTr  (29) 

The  requirement  that  tha  waveform  ands  with  nagativa  alopa  is  mat 
by  restricting  n  to  only  odd  values  and  solving  aquation  (29)  for 
TaLFM  such  that 


T.  .  <2n  f 

c  2 


(30) 


Next,  solve  aquation  (30)  for  n 

n  -  — * 


(o)e  +  2f)  _  i 

2tt  2 


(31) 


For  TaLfM  approximately  aqt«al  to  Td,  Td  is  substituted  into 
aquation  (31)  in  place  of  Ta^p^  and  nLpM  is  chosen  as  the  closest 
integer  (CINT)  to  n  resulting  in 


aKf 

(a)  +  --±)  Td 

C  £ 

2tt 


"UM  -  CIm  1 


(32) 
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Now,  using  in  place  of  n  in  equation  (30)  yields 


(2n_  +  1)tt 

T#LFM  "  sEf - 

™  03  +  (=f) 


(33) 


For  V-chirpod  linear  FH,  Ta  (T«VLFf()  is  found  by  svalusting 
equation  (18)  at  t  ■  T*yLFn  *nd  setting  it  equal  to  1  so  that 
the  waveform  ia  symmetric  about  Ta^p^  and  there  is  no  sudden 
phase  reversal.  No  restriction  is  placed  on  n  in  this  case 
because  T»yLpj|  can  occur  at  any  peak.  In  a  similar  manner  aa 
above,  except  now  Td/2  replaces  Ta  whan  solving  for  n,  nyLpg  can 
be  shown  to  be 


aKf 

(W  +  -"2  )  Td  1 
nVLFM  *  CINT  I— 2 - 4^ - V 


(34) 


and  ^aVLpM  is  found  aa 


t*vlfm 


(2nVLFM  +  2>* 


(U3  + 

c 


-) 


(35) 


The  value  of  Ta  for  the  new  non-linear  FN  waveform  is  found 
when  equation  (16)  ia  normalised,  the  initial  phase  net  to  -it/2, 
and  Ta  substituted  for  t,  yielding 


Xj^CTa)  -  sin  [  (w,  +  0.43  aKf)  Ta] 


(36i 
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Following  tho  same  procedure  used  for  lineer  FM,  the  velues  of  Te 
end  n  for  up  or  down-chirped  non-lineer  FH  can  be  shown  to  be 

(w  +  0.43  aK  )  Td  . 

"SLIM  '  CINT  - U - 1 - P  (37> 


TaVNLFM 


.  (2nVNLFM  +  2^ 

<W  +  0.43  aK_) 
c  t 


In  the  up  or  down-chirped  waveforms,  the  value  n+1 
represents  the  total  number  of  nulls,  (n+l)/2  is  the  total  number 
of  peaks,  and  (n-l)/2  is  the  total  number  of  valleys.  For  the 
V-chirped  waveforms,  4n+2  represents  the  total  number  of  nulls, 
2n+l  is  the  total  number  of  peaks,  and  2n  is  the  total  number  of 
valley? , 

Now,  knowing  that  the  first  critical  point  (a  null)  is  at 
t«0,  the  functior  -is  increasing  in  frequency  and  the  next 
critical  point  .  a  peak,  the  Walkerdid  algorithm  can  be 
implemented. 


41 


Step  It 

Cell  the  location  of  the  currently  known  critical 
point  t^. 

Step  2 t 

Find  the  inatantaneoua  frequency  at  t-  and  call 
it  fx. 

Step  3t 

Calculate  t.  using  t.  •  (t  +  l/4f .) .  Since  the 
frequency  is  increasing*  t.  will  bei»lightly  past 
the  next  critical  point. 

Step  4: 

Find  the  instantaneous  frequency  at  t^  and  call 
it  f 2« 

Step  5: 

Calculate  t.  using  t.  ■  (t.  +  l/4f_).  Since  f 
is  higher  tnan  f.,  t^  will  occur  slightly  before 
the  next  critical  point. 

Step  6: 

Calculate  the  argument  of  the  chirped  waveform  at 
t  and  at  t  and  call  than  arg2  and  arg  , 
respectively.  * 

Step  7: 

Calculate  y  and  y_.  If  the  currently  unknown 
critical  point  is  I  peak  or  valley  use  y.  > 
cos(arg  ),  y  ■  cos(arg  )  otherwise  use  y  * 
sin(arg2),  y^«  sin(arg.7.  The  nulls  of  the 
wave for*  occur  when  the  function  goes  to  xero, 
the  peaks  and  valleys  occur  when  the  derivative 
of  the  function  goes  to  zero. 

Step  8: 

Using  linear  interpolation  calculate  the  value  t 
from  the  equation  t  -  t  -  y,(t,-t  )/(y  -y  ) . 

This  step  essentialy  "homes"  in“on  the  currently 
unknown  critical  point. 

Step  9: 

Calculate  the  argument  of  the  chirped  waveform  at 
t  and  call  it  arg. 

Step  10: 

Calculate  y.  If  the  currently  unknown  critical 
point  is  a  peak  or  valley  use  y  *  cos(arg), 
otherwise  use  y  ■  sin(arg). 

Step  lit 

If  y-0  go  to  step  15. 

Step  12: 

If  the  sign  of  y  equals  the  sign  of  y  ,  let  y  *y 
and  t2«t. 

Step  13: 

If  the  sign  of  y  equals  the  sign  of  y  ,  let  y  »y 
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Step  14 i  Go  to  stop  8. 

Stop  15 i  Store  t  oo  the  current  crlticel  point.  If  t  i> 
leaa  than  To  then  goto  atop  1. 

Step  16 i  If  the  deaired  critical  pointa  were  for 

down-chirped  or  V-chirpud  waveforaa,  uae  the 
ayaaaetry  equationa  (25)  end  (26)  to  adjuat  then. 
End. 


Structural  Layout  of  Tranaducera 


The  trenaducer'a  geonetry  and  the  waveform  are  related  by 


the  acouatic  velocity*  v  ,  of  the  aubatrate.  The  wavelength, 
X,  can  be  calculated  from 


(41) 


The  correapondence  between  the  waveform  end  atructural  layout  of 
the  diaperaive  trenaducer  uaed  in  the  deaign  ia  ahown  in  Figure 
28.  Double  electrodea  were  uaed  to  help  reduce  reflectiona.  The 
x-poaition  of  each  electrode  correaponda  to  the  nidpoint  between 
critical  pointa.  In  order  to  maintain  a  502  duty  cycle 
(approximately  equal  apacing  and  electrode  widtha)  the 
inatantaneoua  frequency,  f  ^ ,  correaponding  to  the  electrode' a 
x-poaition  nuat  be  calculated.  The  widtha  of  the  electrodea  are 
^i/8,  where  ^i  repreaenta  the  inatantaneoua  wavelength. 

In  order  to  ecccount  for  the  frequency  dependence  on  the 
output  of  the  SAW  device,  the  trenaducer  nuat  be  properly 
apodized.  For  an  even  created  waveforn,  this  apodization  ia 
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.V2 

dependent  on  f^(t)  ,  thet  ia  the  amplitude  of  output  la 

inveraely  proportional  to  the  inatantaneoua  frequency  (Her teen, 
Bell  and  Rosenfeld  1973)*  The  largest  aaount  of  overlap,  the 
acouatic  beaawidth  H  ,  occura  at  the  lowest  frequency  and  was 
choaen  ea  100  ^  (the  wavelength  at  the  center  frequency).  The 
aaK>unt  of  overlap  corraaponding  to  each  peak  and  valley  wea 

_V/o 

found  uaing  f.  ,  then  nornalising  to  W#.  These  overlaps 
were  then  used  to  deteraine  the  y-poaition  and  height  of  each 
electrode.  Figure  29  shows  the  effect  of  the  apodisation. 

In  order  to  Maintain  30Z  Metallisation,  dummy  electrodes  are 
used.  These  duamy  electrodes  were  placed  such  that  the  electrode 
gaps  were  ^/g.  Bus  bar  heights  were  choaen  large  enough  (10 
nils)  to  allow  for  bonding  without  having  to  add  additional 
bonding  pads.  Figure  30  shows  the  transducer  with  duauay 
electrodes. 

The  width  of  the  non-dispersive  transducer  was  chosen  to  be 

1.5  X  .  The  distance  between  transducers  was  50  ails  and  a  10 
o 

nil  wide  bar  was  placed  aidway  between  transducers  to  help  reduce 
RF  feedthrough.  Figure  31  shows  the  complete  SAW  device. 

The  FORTRAN  source  code  used  in  the  design  of  the  devices  is 
provided  in  the  Appendix.  The  files  created  are  in  SAWCAD 
STRUCTURE  format  (Nalocha  and  Richie  1984).  The  variable  ISNUM 
is  the  dimension  of  the  X,  Y,  W,  A,  IREP,  XD,  and  YD  arrays  where 

X  «  horisontal  position  of  the  lower  left-hand  corner 
of  the  rectangle 


•  * 

•  • 
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Figure  29.  The  effect  of  apodisatlon. 
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Y  ■  vertical  position  of  the  lover  left-hand  corner 
of  the  rectangle 

W  ■  width  of  the  rectangle 

B  «  height  of  the  rectangle 

A  a  angle  orientation  (always  0) 

IREP  a  nun bar  of  incremental  repetitions  (always  1) 

XD  a  incremental  step  in  the  X  direction  (always  0) 

YD  a  incremental  step  in  the  Y  direction  (always  0) 

The  devices  were  designed  for  a  center  frequency  of  70  MBs, 
e  20  MBs  bandwidth,  a  3  usee  dispersion  time  end  for  fabrication 
on  Y-cut  Z-propogating  lithium  niobate  (501  metallisation 
acoustic  velocity  a  3448  m/s).  The  devices  were  approximately 
12.003  ami  in  width  and  5.483  mm  in  height. 


CHAPTER  VI 


DESIGN  IMPLEMENTATION 

Onc«  th«  STRUCTURE  file  hat  bean  created,  a  photomask  must 
be  produced  ao  that  the  devices  can  be  fabricated.  Though  the 
facility  for  producing  photomasks  is  not  available  at  the 
University  of  Central  Florida,  local  industry  (in  this  case, 
SAWTEK  of  Orlando)  often  provides  the  necesaary  aupport. 

Photomask  Generation 

To  make  full  uae  of  the  photomask,  the  STRUCTURE  files  were 
stepped  and  repeated  until  all  of  the  available  space  was 
utilised.  The  final  product  had  four  copies  of  each  of  the  four 
designs  (linear  FM,  non-linear  FM,  and  the  V-chirped  versions  of 
each)  plus  some  additional  designs  by  other  students.  The  format 
of  the  STRUCTURE  data  had  to  be  converted  to  a  format  known  as 
ELECTROMaSK  in  order  to  be  understood  by  the  pattern  generator, 

ELECTROMASK  data  uses  a  body  centered  coordinate  system, 
assumes  dimensional  data  is  in  units  of  tenths  of  microns  and 
angle  data  ia  in  tenths  of  degrees,  and  has  maximum  aperture 
limitations.  The  subprogram,  SBREAK2,  was  written  to  convert 
general  STRUCTURE  data  (to  include'  any  arbitrary  angle  and  any 
number  of  repeats)  to  ELECTROMASK  data  and  is  included  in  the 
Appendix . 
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Fabrication 

Tha  fabrication  of  davicas  was  accomplished  using  the  new 
clean  room  facilities  at  the  University  of  Central  Florida.  For 
a  detailed  explanation  of  the  fabrication  process  see  Vigil  and 
Yapp  (1984).  Briefly,  the  process  consists  of  the  following 
steps s 

1)  Clean  the  lithium  niobate  wafer. 

2)  Deposit  aluminum  using  flash  evaporation  or  sputtering 
techniques . 

3)  Apply  photoresist  to  the  metallized  surface  and  soft 
bake. 

4)  Use  the  mask  aligner  to  make  contact  between  the  wafer 
and  mask,  then  expose  to  the  ultraviolet  source. 

5)  Develop  the  wafer  to  remove  the  exposed  photoresist. 

6)  Place  the  wafer  in  an  aluminum  etch  to  remove  the  desired 
metal . 

7)  Dice  the  wafer  to  separate  the  individual  devices. 

8)  Mount  the  devices  to  a  header. 

9)  Bond  the  devices  to  the  header. 

10)  Apply  absorbing  material  to  the  ends  of  the  devices  to 
reduce  edge  reflections . 

A  few  problems  were  encountered  during  fabrication  which 
would  affect  the  performance  of  the  devices.  One  such  problem 
was  the  inability  to  achieve  a  smooth  layer  of  photoresist  in 
step  3.  This  resulted  in  a  number  of  fingers  on  the  devices 
being  open  or  shorted.  Another  problem  was  created  back  when 


the  STRUCTURE  data  was  stepped  and  repeated.  There  should  have 
been  enough  distance  left  between  the  devices  to  allow  for  a 
sufficient  amount  of  absorbing  material  to  be  placed  at  the  ends 


of  the  device 


CHAPTER  VII 


CONCLUSIONS 

Fundamental  radar  concepts  have  been  reviewed,  emphasizing 
the  need  for  matched  filtering  and  the  desirable  characteristics 
of  the  transmitted  pulse  and  its  autocorrelation  function.  A  new 
non-linear  FN  pulse  compression  technique  has  been  reviewed  using 
linear  FN  as  the  basis  of  comparison,  A  simulation  of  the 
Doppler  effect  on  the  matched  filter  output  has  beer, 
accomplished,  including  the  effects  on  V-chirped  signals.  SAW 
device  design  considerations  have  been  presented  and  an  efficient 
algorithm  for  finding  the  critical  points  of  a  chirped  waveform 
was  introduced.  Computer  software  was  written  and  is  presented 
in  the  Appendix  to  implement  the  SAW  device  design  and  the 
subsequent  conversion  of  data  for  photomask  generation.  SAW 
device  fabrication  was  briefly  reviewed. 

The  Doppler  simulation  showed  that  the  new  non-linear  FM, 
unlike  linear  FM,  could  be  used  as  a  moving  target  indicator  by 
examining  only  the  matched  filter  output  in  the  time  domain. 
Furthermore,  the  new  non-linear  FM  offers  more  dynamic  range  than 
the  linear  FN  because  the  time  shift  is  less  than  that  of  linear 
FM. 

In  conclusion,  it  is  recommended  a  new  photomask  be  produced 
which  provides  a  layout  allowing  for  a  suf  *‘.cient  amount  of 


absorbing  material  to  ba  plaead  at  tha  ends  of  tha  davicaa. 

Also,  a  simulation  which  includes  the  combination  of  tha  naw 
non-linear  FM  with  windowing  and  tha  Doppler  affect  may  add  to 
tha  resolution  capabilities  of  tha  pulse.  Further  research  could 
include  tha  naw  non-linaar  FM's  potential  in  a  multiple  target 


environment . 


APPENDIX 


COMPUTER  PROGRAMS 


o  u  if  uo 


prograa  chirpy 


data  of  last  revision: 


09-16-86 


for  aore  information  contact: 


O.C.  Ma lochs 
J.C.  Walker 


c 


prograa  chirpy 
character*!  ask 

coaaon/patl/  x(6000),y(6000),w!6000),h(6000),a<6000), 
2  inp(6000),xd(6000),yd(6000),isiua,lref 


c 

500 


600 


700 


**  HEMJ  ** 

call  cla 
write!*,*)  • 
write!*,*)  ' 
write!*,*)  ' 
write!*,*)  • 
write!*.*)  ' 
write!*,*)  ' 
write!*,*)  ' 
write!*,*)  ' 
write!*,*)  ' 
write!*,*)  ' 
write!*,*)  ' 
write!*,*)  ' 
write!*,*)  ' 
write!*,*)  ' 
write!*,*)  ' 
write!*,*)  1 
write!*,*)  ' 
write!*,*)  ' 
read!*, 600)  ask 
foraat(al) 

if  (ask.eq. 'c'.or.ask.eq. 
if  (ask.eq. 'p' .or. ask.eq 


Dispersive  S/M  Device  Design' 


(C]reate  File' 


(Pjrint  File' 
(S)ave  File' 
iQlult' 


Enter  Choice' 


if  (ask.eq  ' 
if  (ask.eq 
goto  500 

stop 

end 


.or.asx.eq 


*C* )  call  create 
'P')  call  pmt 
S')  call  savlt 


or. ask.eq. 'Q' )  goto  700 


c 


subroutine  els 

do  800  1*1,22 
write!*,*)'  ' 
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800  continue 
return 


subroutine  create 
character*!  aek 

coaeon/patl/  x(6000),y(6000),v(6000),h(6000),a<6000), 
2  lrep(6000)txd(6000),yd(6000),lsnue,lref 
coneon/cheq/  akf ,arg,fl,lfa,pi,t,ta,wl,wl 
real  tlae(3000),ov<3000),a!6000> 

c  **  MW  t  VERIFICATION  ** 

1000  call  da 

1100  write!*,*)  *  ' 

write!*,*)  '  1)  Linear  FM' 

write!*,*)  i)  Non-Unear  FM' 
write!*,*)  '  ' 

write!*,*)  1  Enter  Choice  (1,2)' 
write!*,*)  '  ' 
read!*,*)  lfe 

lf(lfe.ne.l.and.lfe.ne.2)then 

write!*,*)  •  Input  error  -  try  again' 

goto  1100 

endif 

1200  write!*,*)  '  ' 

write!*,*)  *1)  Up  or  Dow  Chirp' 
write!*,*)  *2)  'V  Chirp' 
write!*,*)  '  ' 

write!*,*)  '  Enter  Choice  (1,2)' 
write!*,*)  '  ' 
read!*,*)  lahape 
write!*,*)  '  ' 

If ( lahape . ne . 1 . and . lahape . ne . 2 ) then 

write!*,*)  '  Input  error  -  try  again' 

goto  1200 

endif 

write!*,*)  '  ' 

write!*,*)  '  Hlnleua  Frequency  (  MHa  )?' 
write!*,*)  '  ' 
read!*,*)  fit 

write!*,*)  '  ' 

write!*,*)  '  Haxlaua  Frequency  (  MHz  )?' 
write!*,*)  '  ' 
read!*,*)  fht 


write!*,*) 


•  I 
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1300 


1400 


write!*,*)  '  Dispersion  T1m  (  umc  )?' 
write!*,*)  '  ' 
mdl*,*)  tdt 


•  • 

Data  in  units  of 


*  i 


1)  wavelengths' 

2)  ailllaeters' 

3)  ■lcroMtsrs' 

4)  alls1 


Enter  Choice  <1,2, 3, 4)' 

it 


vrite!*,*) 
write!*,*) 
write!*,*) 
write!*,*) 
write!*,*) 
write!*,*) 
write!*,*) 
write!*,*) 
read!*,*)  lscale 

if ( 1 scale. It . 1 .or . I  scale,  gt .4 ) then 

write!*,*)  '  input,  error  -  try  again* 

goto  1300 

endlf 

write!*,*)  '  ' 

write!*,*)  1  SAW  Velocity  (  eeters/aec  )?* 
write!*,*)  '  ' 
read!*,*)  va 


Up  or  Down  Chirped' 
"V"  Chirped' 


call  els 

lfdahape.eq.Dwrlte!*,*) 

If ( 1 shape . eq . 2 )  write ( *,* ) 
write!*,*)  '  ' 
lf(ife.eq.l)write(*,*)  '  Linear  FN' 
if (ifa.eq.2)wrlte(*,*)  '  Non-Linear  FM' 
write!*,*)  '  ' 

. Mlnlaus  Frequency  -  *,flt,'  We' 


Maximal  Frequency  -  ' , fht , '  MHz' 
Dispersion  Tine  -  ' ,tdt, '  usee' 


write!*,*) 
write!*,*) 
write!*,*) 
write!*,*) 
write!*,*) 
write!*,*) 

if  liscale.eq.l)  write!*,*) 
if  (iacale.eq.2)  write!*,*) 
if  (lscale.cq.3)  write!*,*) 

If  !iscale.eq.4)  write!*,*) 
write!*,*)  '  ' 
write!*,*)  '  SAW  Velocity  ■ 
write!*,*)  '  ' 
write!*,*)  '  * 

write!*,*)  *  Is  this  correct  (y/n)?* 
read!*, 1400)  ask 
foraat(al) 


'  Units  -  wavelengths' 
'  Units  -  alllineters' 
'  Units  •  alcroaeters' 
'  Units  •  alls' 

',va,'  aeters/sec' 


if  (ask.eq. 
write!*,*) 
write!*,*) 
write!*,*) 


'n'.or.ask.eq.'N')  goto  1000 


PROCESSING" 


c  **  INITIALIZATION  ** 


fi-4.0*atan(1.0) 
L-flt*1.0e6 
fh-fht*1.0e6 
td-tdt*1.0*-6 
if  (lahapa.aq.2)  td-td/2.0 
akf-2.0*pl*(fh“fl) 
wl-2.0*»i*fl 
wh-2.0*pi*fh 


If  (Ifsuaq.l.and.lshapa.eq.l)  than 
n-td*(0.3*fl*0.5*fh) 
ta-(2.0*ml.0)/(fl*fh) 
andlf 

if  (lfa.eq.l.and.lshaca.eq.2)  than 
n-td*(0.5*fl«0.5*fh)»0.25 
ta-(2.0*n+0.5)/(fl+fh) 
andlf 

if  (lfa.*Q.2.and.iahap*.eq.l)  than 
n-td*?0.57*fl+0.«*fh) 
ta-(2.0*h+1.0)/(1.14*fl*0.86*fh) 
andlf 

if  (lfa.cq.2.and.iahapa.cq.2)  then 
n-td*(0. 57*f  l*0.43*fh)  *0.25 
ta-(2.0*n+0.5)/(1.14*fl+0.86*fh) 
andlf 


If  (ishape.eq.l)  then 
ncp-4*n+3 
nf-4*n+2 
andlf 

If  (lshape.eq.2)  then 
ncp-4*n+2 
nf-«*n*2 
andlf 


lsrua«2*nf*17 
fo-(fl*fh)/2.0 
rlo-va/fo 
height-100. 0 
wa -height *rlo 
bbh-0.000254/rlo 
d-0.00127/rlo 
utw-1.5 


!  Acoustic  Beaawidth  in  wavelengths 

!  bus  bar  height  is  10  ails 

!  50  ails  between  transducers 

!  unweighted  transducer  is  1.5  wavelengths  wide 


**  CALCULATE  NULLS,  PEAKS  &  VALLEYS  ** 

tla*(l)-0.0 

tl-0.0 

do  1700  1-2, ncp 
r-l 
t-tl 

call  fraq 
fl-fi 

t2-tl+17.0/64,0/£l 


i 


* 
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t-t2 

call  freq 
f2-fl 

call  arguaant 

It’  Ur/5.0 )  .ft. (1/2))  ypos»aln(arg) 

It  ((r/2.0).eq.(l/2>)  ypoe-coa(arg) 

t»t3 

call  argument 

It  ((r/z.0).gt.(l/2>)  ypoa>aln(arg) 

If  ((r/2.0).eq.(l/2)i  ypoe-coe(arg) 

).*q.*gn(ypo$3))  than 

wrltel*,*)  'ERROR:  design  apaclflcatlona  cannot  be  wet!' 

try  althar  a  longer  dispersion  tlae' 
or  sataller  bandwidth.' 


1499 

1500 


1550 

1600 

1700 


1800 


write!*,*) 
write!*,*) 
write!*,*) 
write!*.*) 
goto  1499 
endlf 

t-t2-ypos2*< t3-t2 ) / (ypos3-ypo*2 ) 

If  (t.eq.t2.or.t.eq.tJ)  goto  1550 
call  araiaant 

If  (!r/2.0).gt.!l/2))  ypos-sln(arg) 
If  !!r/2.0).eq.!l/2l)  ypowcoeiarg) 
If  (sgn(ypoa).eq.O)  tl-t 
If  (s|n(ypos).eq.sgn(ypos2))  then 

ypos2-ypos 

endlf 

If  (s|r\typos).eq.sgn(ypos3))  then 

ypoa3-ypoa 

endlf 

If  (vpos)  1500,1600,1500 
tl-t2 

If  (abs(ypos3).lt.abs(ypos2) )  tl«t3 


ENTER  <  CNTRL  >  Y  to  return  to  systee' 


tlne(l)«tl 
continue 

If  (lshape.eq.1)  goto  1900 
do  1800  l»l,(ncp-l) 

tt*e!ncp*l)»2.0*ta-tlee(ncp-i) 

continue 


use  syewetry  for  ’V  Chirp 


c  **  WERLAPS  ** 

1900  oveax*0.0 

do  2000  1*2, nf ,2 
t-tlne(i) 

If  (lshape.eq.2.and.t.gt.ta)  t*2.0*ta-t 

call  freq 

ov( 1/2 )-f !**<-!. 5) 


o  u  u  o 


2000 


if  <Ov(l/2).tC«(MMW)  «mmx-ov(1/2) 

continue 


2100 


continue 


C  ■**  WIDTHS  *  x  positioning  ** 

do  2200  1-l.nf 
tl-t 

if  (lahape.eq.2.and.e.gt.ta)  t-2.0*ta-t 
call  fraq 

»<l)-1.0/B.0/fi*fo 

w<nfH)-w(l) 

x(i)-tl*fo-w(l>/2.0 

x(nf*l)«x(l) 

2200  continue 


C  **  HEIGHTS  *  ¥  POSITIONING  ** 


2250 


2300 


...  .  ^  !  irk  is  the  distance  between  hue  bora 

h(lMrk*otf(l)*bfah)/2.0 
h(nM)-rk*bbh-h<l)-w{l) 
y(l)-rkO.O*t*h/2.0-h(l) 
y(nfn)-bbh/2.0 
do  2300  l-2,(nf-2>.2 

h(l)«rk*ov(i/2)+bbh~h(i-l) 

h(nf+l ) -rk*bbh-h( l)-w(i) 

h(l+l)-h(l) 

h(nf *1+1 )*h(nf  ■ *1 ) 

r-1 

if  <<r/4.0),et.(l/4))  than 
y(i)*bhh/2.0 

y(nf+i)«rk+3.0*bbh/2.0-h(nf*l) 

endif 

if  <<r/4.0),eq.(i/4))  then 
y(l)-rkO.O*bbh/2.0-h(l) 
y(nf«l)»bbh/2.0 
endif 

y(l*l)-y(l) 

y(nf*l*l)»y(nf+l) 

continue 

h(nf )>rk«ov{nf/2)+bbh-h(nf-l) 

h(2*nf  )-rk*bbh-h(nf  )-w(nf) 

y(nf)>bhh/2.0 

y(2*nf )-rk+3.0*bfch/2.0-h(2*nf ) 


Make  sure  it’a  ok.  If  an  olectrodo  leaa  than  half 
«*  *■  bar  height  than  we  need  oore  distance  between 
the  bus  bars  to  nake  it  lnpleeentable  (though  this 
does  not  guarantee  a  practical  design). 
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MOO 


lfUe-0 

do  MOO  l«l.!2*nf) 

If  (h(l).lt. (bbh/2.0))  than 

lfU|-lflarl 

continue 

if  (lflaf.at.O)  than 
rk-r*»1.0 

If  (Tk.gt.f.Otwlaht)  then 

write!*,*)  'distance  between  bue  bare  la  now:' 
write!*,*)  rfc.  '  wavelengths' 
wrttet*,*)  '  1 

write!*,*)  'acouatlc  be —width  la:1 
write!*,*)  hel^it,  'wavelengths' 
write!*,*)  •  • 

write!*,*)  'Enter  <  CHIXL  Y  >  If  you  went  to  abort* 
endlf 
goto  2230 
endlf 


c  **  BUS  BARS  ** 

x(2*nf*l)»  -0.23 
x(2*nf+2)-  -0.25 

If  (lahapa.eq.l)  w!2*nf»l)»ta*fb*0.5 

If  (lahape.eq.2)  w!2hnf«l)«2.0*ta*fo*0.5 

w(2*nf*2)-w(i*nf*l) 

h(2*nf*l)-bbh 

h(2*nf*2)-bbh 

y(2*nf*l)-0 

y(2*nf+2)-rk*bbh 


C  **  NCN-01SKKSIVE  TRANSDUCER  ** 

x!2*nf*3>-v(2*nfd)*d*0.8 
x(2*nf44)-w!2*nt*2)*d*0.8 
w!2*nf*3)-utw*d*0.4 
w!2*nf«4)«utw*d*0.4 
h!2*nf*3)-bbh 
h(2*nf*4)-bbh 
y(2*nf*3)-0.0 
y(2*nf+4)«rk*bbh 
do  2300  l-l,(utw*4) 
rl-1 

x(2*nf*4d)»x(2*nf'*3)'*l.Q/16.0*(rl-1.0)/4.0*d*0.2 
w( 2*nf +4*1 )«1. 0/8.0 
h(2*nf  ♦4*l)«(rV*helaht*bbh)/2.0 
If  (l.eq.l.or.l.eq.4.or.l.eq.3)  than 

y(2*nf*4d)-rk*3.0*bbh/2.0-h(2*n£*4+l) 

endlf 

If  (l.eq.2.or.l.eq.3.or.l.eq,6)  than 
y(2*nf*4*l)-bbh/2.0 


andlt 

x(2*nf»10»i)-xt2«nf«(*l> 

%H2%«f  ♦10*0-1. Q/t.O 
h(2*nf *10*1  )-rtt*fettv-h(2*nf*4«l  )-l  .0/8.0 
if  U.M.l.or.l.eq.d.or.l.M.S)  then 
y<2*nf*l0*l)-tth/2.0 
endlf 

it  (l.«.i.or.i.iq.).or.t.M.(l  than 

y(2(nf *10*1  )-rtt*3.0»wih/2.0-h(2*itf  *10*1 ) 
•ndif 

2900  oontlnut 

x(im»-2)>M(M«l)40.i*d 
wdenw2)-d/5.0 
h( ianun-2 )-rit*2.0*fcbh 
y(lenue-2)-0.0 


c  **  SOTTING  ** 

do  2(00  l-l,(lenun-2) 

•(i)-x(l) 

2(00  continue 

x(l)-e(2*nf*l) 

x(2>-e(2*nf*3> 

x(3)-«<2*n£*4> 

x<4)-e(2*nf*2) 

io  2700  1-5  ,(2*nf*3),2 

2700  continue 

JTmoo  l-(,<2*nf*4),2 

iSKe(j) 

2800  continue 

x(2*nf+5)»etlenun-2) 

J-2*nf*A 

do  2900  l-(2*n£+(),(ienu*-3),2 

x<ft-e(j) 

2900  continue 

J-2*nf«4*Utw*4.0 

do  3000  l-(2*nf*7),<l*"»-2),2 

3000  continue 

do  3100  l-l,(len»-2) 

•(l)-y(t) 

3100  continue 

y(l)-e(2*nf«l) 

y(2).»(2*nf*3> 

y(3)-e(2*nf*> 

y(4)-e(2*nf*2) 

>0 


do  3200  l-5,(2*nf+3),2 
J-J*l 
y(I)-s(j) 
continue 
J-nf 

do  3300  l«6,(2*nf+4),2 
continue 

y(2*nf+5)*s(lsnun-2) 

1-2*nf+4 

oo  3400  i«(2*nf+6),(lsnu»-3),2 
J-J*l 
y(i)-«(j) 
continue 
J-2*nf+4+utw*4.0 
5°  3500  l“(2*nf+7),(lnmj»-2),2 
J-j*l 
yU)-8(J) 
continue 

do  3600  1-1, (isnum-2) 

*(l)-w(i) 

continue 

w(l)-s(2*n£+l) 

w(2)-s(2*nf+3) 

w(3)-s(2*nf+4) 

w(4)-s(2*nf+2) 

J-0 

do  3700  l“5,(2*nf+3),2 
>j+l 
w(l)-s(M 
continue 
j-nf 

do  3800  l-6,(2*nf+4),2 

w(i)-s( J) 
continue 

w( 2*n£+5 )-s ( lsnun-2 ) 

J-2*nf+4 

ao  3900  l-(2*nf+6),(lsnum-3),2 

J-J+l 

w(l)-s(J) 
continue 
J-2*nf+4+uLw*4 . 0 
do  4000  i-(2*nf+7),(lsnun-2),2 

w(i)-s(j) 

continue 

do  4100  1-1,1 lsnum-2) 
sfi)-h(l) 
continue 
h(l)-8(2*nf+l) 
h(2)-s(2*nf+3) 
h(3)-s(2*nf+4) 
h(4)-sf2*nf42) 


4200 

4300 

4400 

4500 


1-0 

do  4200  l-5,(2*nfO),2 

fSlSa(J) 

continue 

J-nf 

do  4300  i«6,(2*nf+4),2 
continue 

h!2*nf+5)-a(isnu»-2) 

J-2*nf-*4 

do  4400  i-(2*nf*6),(ienu»-3),2 

i!<  &•<.)) 

continue 

£2*nf+4mtw*4.0 

4500  i»(2*n£+7),(lanun-2),2 

iM)-a(j) 

continue 


c  **  SCALING  ** 

if  (lecale.cq.l)  ef-1.0 
if  (iscele.eq.2)  af-rlo*1000.0 
if  (iscale.eq.3)  af-rlo*1.0e6 
if  (lacale.eq.4)  a£-rlo*1.0e5/2.54 
do  4600  i-l,isnua 

x(l)-(x(l)-xnin)*ef 
y(  l  My  ( i )-ynin)*af 
w(l)^(l)*af 
h(i)-h(l)*a£ 

4600  continue 

4700  do  4800  i-l.ianun 

•(i)-0.0 
irep(l)-l 
xd(  U-0.0 
yd(i)-0.0 
4800  continue 

return 

end 


aubroutine  freq 

Ihia  aubroutine  finds  the  inatantaneoua 
frequency  of  the  chirped  waveform 

coMaon/cheq/  akf , erg, £1, ifn.pl, t,te,wi,wl 


c 

c 
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if  (ifa.sq.l)  wi-wi+»kf*t/t« 
if  (lfm.cq.2)  then 

wi-vl+akf *(0.43+0. 5*eoe ( pi*( t-te )/ta ) 
2  ♦0.07*eoe(2.0*pi*<t-ta)/te>) 

endlf 

fi-wi/2.0/pl 

return 

end 


subroutine  arguaent 

Subroutine  calculates  the  argument  of  chirp  waveform 
equation.  USed  by  calling  prograe  to  find  critical 
points.  If  seekity  a  null  of  the  waveform,  sln(arg) 
is  used.  If  seeking  a  peak  or  valley,  the  derivative 
of  the  function,  cos(arg),  is  used.  Either  way,  the 
calling  program  la  looking  for  the  overall  result  to 
go  to  aero. 

coaaaon/cheq/  akf ,  arg,  f  1 ,  lfa  5pi ,  t ,  ta  ,wi  ,wl 

if  (lfn.eq.l)  then 

erg-wl*t*akf/2 . 0/ta*t*t 
endlf 

if  (lfs.eq.2)  then 

arg«*»l*t*akf*(0.43*t+0. 5*ta/pi*sin(pl*<  t-ta  )/ta  )♦ 
2  (0.07*ta)/(2.0*pi)*sln(2*pi*(t-ta)/ta)) 

endlf 

return 

end 


integer  function  sgn(x) 

if  (x.eq.0.0)  sgn-0 
if  (x.gt.0.0)  sgn-1 
if  (x. It .0.0)  sgn«-l 

return 

end 


subroutine  pmt 

coamnn/patl/  x(6000),y(6000),w(6000),h(6CX)0),-(6000), 
2  lrep(6000),xd(6000),yd(6000),lsnun,lref 


return 


* 


c  this  is  just  s  Modified  version  of  subroutine  vriteo 

subroutine  sevit 

cooeon/patl/  x(6000),y(6000) ,w(6Q00) ,h(6000) ,a<6000) , 
2  lrep(6000),xd(6000),yd(6000>,isnus,iref 

character  filout*10 


c  terminal  input 

100  write<6,*)  '  «<  SYSTEM  WRITE:  »>• 

write(6,*)'  ' 
wrlte(6,1005) 

1005  fornat(x, 'ENTER  output  file  nane  :  — >  '  ,$) 

read (5, 1000, err- 10)  filout 

1000  foreat(alO) 

open( 10,  flle-filout, status- 'unknown' ,err-10) 
close( 10, status- 'delete' ) 
open(10,file-fllout, status- 'unknown' ,err-10) 
goto  20 

10  write(6,*)'  ***  ERROR  in  file  nane  -  try  again  ***' 

goto  100 


c  write  out  file  data 


wrlte(l0,2000)  icont 

2000  foreat(x, 'icont-' ,15) 

write< 10,2001)  itype 

2001  fornat(x,'ltype-',lA) 


write (10, 2002)  fo 
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2002  format (x,'fo  -\el7.9) 

wrlte!10,2003)  tflo 

2003  format!x,'tflo  -\el7.9) 

write! 10, 2006)  tfhl 
2006  format (x, 'tfhl  -',el7.9) 

write! 10, 2005)  num 

2005  format(x,'num  -',14) 

write! 10,*)  aero, aero 
write!  10,*)  aero,  sera 

write!  10,2006)  Unue 

2006  format  (x, 1  lenum-'  ,15) 

write! 1C, 2007)  lref 

2007  fonaattx, ‘lref  -',11) 

do  400  l-l,i»nue 

write! 10,*)  x(l),y(l),w(l),h(l),a!l),)rep(l),xd(l),yd(l) 
400  continue 


cloee(lO) 

return 

end 


c 


ouu 
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subroutine  sbreak2( raw rk,n£ lash) 


c  - - - - - - 

For  Information  contact:  D.C.  Malocha 
S.H.  Richie 
J.C.  Walker 

c  - - - 

c  Date  of  last  revision:  Aupat  19,  1986 

c  - - - - - - 

c  This  subroutine  is  used  when  generating  Magtapes  using  the 
c  ELECTROMASK  (SA WIHC)  format.  Structure  data  (body  centered,  mm) 
c  is  taken  and  new  rectangles  are  generated  to  eliminate  repet  1- 
c  tlons.  Angle  data  is  resolved  such  that  it  is  in  **-»  range  from 
c  0  to  90  degrees;  swapping  height  and  width  data  if  necessary, 
c  Rectangles  are  broken  up  Into  smaller  rectangles  if  necessary 
c  to  satisfy  maximum  aperture  limitations. 


integer  nrec(99) 
character  chstr*ll 
real*8  chnun 
byte  sbuf( 512,2000) 

common/ chconv/  chstr  .chnum 

common/ patl/  x06000)1y(36000),w(36000),h(36000),a(36000), 
2  lrep( 36000) ,xd< 36000), yd< 36000), Isnum.lref 

common/buffer/  sbuf 

comm on/add reas/  numrec,numbyte,nrec,ntape 
common/flash/  xg,yg,wg,hg,ag,lrg,xl,xain,xmax,ymln,ymax 


c  maxima  aperture  for  SAWTEX  Is  1.524  am 

apmax-15240  !  units  are  tenths  of  micrometers 

pi”4*atan(1.0) 


c  load  first  begining  of  record  symbol 

if (numrec.gt .l.or.numbyte.gt .l)then 
numbyte-1 
numrec-numrec+1 
end  If 
lstrlen-1 


*- 


67 


* 

4 


ch*tr(l:l)-'<* 
call  ldabuf(lstrlen) 
neark-n*ark+l 

c  anal*  data  la  rcaolvad  by  firat  putting  the  data  in  the  range 

c  fro*  -180  to  180  degreea  then  adjusting  the  angle  (and  height 

c  and  width  data  if  necessary)  such  that  all  angle  data  is  in  the 

c  range  froai  0  to  90  degrees  (angle  data  is  in  tenths  of  degrees) 

do  1000  1-l.isnua 

300  lf(a(l).gt.l800)a(l)-a(i)-3600 

lf(a(i).lt.-1800)a(l)*a(i)+36G0 
lf(a(l).gt.l800.or.a(i).lt.-1800)goto  500 
if(a(l).It.0)a(i)-a(iM800 
if(a(i).ge.900)then 
a(i)-a(i)-900 
ww-w(l) 
w(i)4i(i) 
h(i)-ww 
endif 
1000  continue 

c  initialise 

do  3000  i-l,lsnue 
wg-w(l) 
hg-h(i) 
ag-a(l) 

agr-a(l)*pi/l800  !  convert  tenths  of  degrees  to  radians 

lrg-lrep(i) 

nwrep-1 

nhrep-1 

wov-0 

hov«0 

lf(wg.gt.apeax)  than 
nwrep-wg/apeax*l 

1100  vov-(rwrep*apeax-vg)/ (nwrep-1) 

if  (wov. It . 10.0 ) then 
nwrep-nwrep+1 

ro  1100 

wg-apeax 

endif 

lf(hg.gt.apeax)then 

nhrep-hg/apeax+1 

1200  hov- ( rhrep*apeax-hg ) / ( nhreo-l ) 

if (hov.lt. 10.0)  then 
rhrep-rhrep+1 

rol2C0 

hg-apeax 

endif 


do  4000  J-O.(irg-l) 


!  repeats 


1500 


do  3000  k»0,(nwrcp-l) 
do  2000  l-0,(t*rep-l) 
xl>x(l)«Md(l)*j 
yt-y(l)*yd(l)*j 
lf(agr.eq.O)then 

x*»xi-w(  i  )/2*wg/2*(wR-vw  )*k 

yt-yl-h(l)/2^/2+(hr-how)*l 

rlSOO 

xlp-xi*co*(sgr )+yi*sln( egr ) 
ylp-yl*cos  ( ur )  -xl*sin(  a*r ) 
xgp-xip-w(  1 )/  2+wg/2+ (wg-wov )% 
ygp*yip-h(  1  )/2+hg/2*  (hg-hov  )*1 
xg-xgp*coa(agr )-ygp*sln(agr ) 
yg-ygp*cos(agr)»xgp*sln(agr) 
lf(nuMbwte.gt.472)then 

abut (nuibyte.nuerec)-’ . • 
•buf  ( 1  ,nuerec+l )» 1  <  * 
nusbyte-2 
r*jir*c-re«rec*l 
endlf 


call  flash2(nsarktnflaah) 
2000  continue 

3000  continue 

4000  continue 
5000  continue 


!  widths 
!  heights 


c  load  end  o£  record  syebol  for  last  rec  of  file 

lstrlen-1 
chstr(l:l)-'.' 
call  ldsbuf(latrlen) 
neark>nurk+l 


return 

end 
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